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While most stars are believed to form in stellar clusters, the formation and early evolution of young stellar clusters is still largely 
unknown. Improving our knowledge of the earliest phases of clustered star formation is crucial for understanding the origin of the 
stellar initial mass function and the efficiency of the star formation process, which both play a key role in the evolution of galaxies. 
Here, we present an analysis of the Aquila Rift complex which addresses the questions of the star formation rate (SFR), star formation 
efficiency (SFE) and typical lifetime of the Class protostellar phase in two nearby cluster-forming clumps: the Serpens South and 
W40 protoclusters. We carried out a 1.2 mm dust continuum mapping of the Aquila Rift complex with the MAMBO bolometer 
array on the IRAM 30m telescope. Using a multi-scale source extraction method, we perform a systematic source extraction in our 
millimeter continuum map. Based on complementary data from the Herschel Gould Belt survey and Spitzer maps, we characterize the 
spectral energy distributions (SEDs) of the 77 millimeter continuum sources detected with MAMBO and estimate their evolutionary 
stages. Taking advantage of the comprehensive dataset available for the Serpens South region, spanning wavelengths from 2 jim to 
1.2 mm, we estimate the numbers of young stellar objects (YSOs) at different evolutionary stages and find a ratio of Class to Class I 
protostars N(Q)/N(I) =0.19-0.27. This low ratio supports a scenario of relatively fast accretion at the beginning of the protostellar 
phase, and leads to a Class lifetime of ~ 4-9 xlO 4 yr. We also show that both the Serpens South and W40 protoclusters are 
characterized by large fractions of protostars and high SFRs ~20-50 M Myr~'pc~ 2 , in agreement with the idea that these two nearby 
clumps are active sites of clustered star formation currently undergoing bursts of star formation, and have the potential ability to 
form bound star clusters. While the formation of these two protoclusters is likely to have been initiated in a very different manner, 
the resulting protostellar populations are observed to be very similar. This suggests that after the onset of gravitational collapse, the 
detailed manner in which the collapse has been initiated does not affect much the ability of a clump to form stars. 

Key words. Stars: formation, circumstellar matter - ISM: clouds, structure - Individual objects: Aquila Rift, Serpens South, W40 
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1. Introduction 

The formation and early evolution of protoclusters 

While most of the young pre-main sequence (PMS) stars in 

molecular clouds are found in groups and clusters rather than 



being uniformly distributed throughout the clou ds (|Testi et al. 
|1999l |Carpenter|2000||Lada & Lada|2003l see also|Bressert et al. 
for a discussion of clustered star form 



2010] for a discussion of clustered star formation), the forma 
tion and early evolution of stellar clusters are still poorly known. 
Improving our knowledge of the earliest phases of clustered star 
formation holds the key to understanding the distribution of stel- 
lar masses at birth since the similarity of the stellar initial mass 
function (IMF) of young stellar clusters with the field stars IMF 
may imply it is linked with the cluster formation process itself 
(|Palla & Stahler|[T999l |Hillenbrand & Carpenter|[2000) |Meyer 
|et al.|2000| >. " 

Protoclusters are parsec-scale gas-dominated cluster- 
forming clumps, which contain a large amount of dense gas 
material (several Jeans masses), and are actively forming 
protostars, as testified by the significant population of prestellar 
cores and embedded Class objects detected toward them 
at (sub-)millimeter wavelengths (see Mott e et al.||1998) |Hurt| 



Based on observations carried out with the IRAM 30m telescope. 



IRAM is support ed by IN SU/CNRS (Fran ce), MPG (Germany), and 
IGN (Spain). 



& Barsony|[T996l |Myers|[200T] |Testi & Sargent] [1998) . While 
protoclusters are the direct progenitors of infrared embedded 
clusters, only a small fraction of them eventually result in bound 
star clusters, as it has been suggested (Lada & Lada|2003j) that 



most embedded clusters are not able to remain bound when they 
rapidly lose their gas mass and evolve into revealed clusters. 
Their stellar component therefore disperse to become field 
stars. This high infant mortality has been proposed to be due 
to insufficient star formation rates (SFR) to fight against fast 
gas loss (Goodwin 2009; Baumgardt & Kroupa 2007), but may 
also rely on the dynamical state of the young stars and left-over 
material during the later evolution of the embedded cluster 
dSmith et al.|20TTT >). 

With ages of only a few ~10 5 yr in some cases (Peretto et al. 



20071, protoclusters are therefore the youngest precursors of 



star clusters. Their youth makes them ideal targets to investigate 
the formation and early evolution mechanisms of clustered star 
formation, since their properties are likely to still reflect the 
physics of their formation. Deep (sub-)millimeter continuum 
observations of cluster-forming clumps can set strong con- 
straints on the early properties and conditions for clustered star 

[19981 |Testi & Sargent[[T998] l, the 



formation (cf. Motte et al 



efficiency of the conversion of gas into stars from well-defined 
mass reservoirs (cf. |Enoch et al.||2007) |Belloche et al.||2011) , 
and the typical lifetimes of the youngest embedded objects, i.e., 
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Class pro tostars ( | Andre et al.|2000| |Visser et al.|2002| |Enoch 



et al.| 2009 >. Moreover, only (sub-)millimeter line observations 
can probe the dynamics of cluster-forming clumps ( |Peretto| 
eTaLl[20fJ6l |Andr6 et al .120071 |Kirk et al.|2007l >, as well as the 
role of feedback from embedded YSOs on the early evolution 



of young protoclusters i 


Knee & Sandell 


2000 Ginsburg et al. 


2009( |Mauryet al.|2009 


Nakamura et al. 


[2011aJ. 



Several observational studies targeting both low-mass and 
intermediate-mass protoclusters have shown that embedded 
clusters form stars more efficiently than global molecular clouds. 
Indeed, typical star formation efficiencies (SFEs) of 10^-0% 
HWilking & Lada|1983[|OTmi & Testi|2002l|Lada & Lada|2003| l 
have been estimated for embedded clusters, while molecular 
clouds have SFEs of only a few percent (Myers et al. 1986 



Evans et al. 2009 ; L ada et al.|2Q10) >. The SFE and star formation 
rate (SFR) are known to vary from region to region, but few stud- 
ies have focused on comparing these properties for protoclusters 
formed from different initial conditions. Moreover, past estima- 
tions of SFEs and SFRs have generally relied on only partial 
population counts and limited source classification, as it is very 
difficult to detect and characterize complete protostellar popu- 
lations with limited wavelength coverage. Indeed, in order to 
probe complete YSO samples from Class protostars to Class II 
objects and properly determine their evolutionary stages, deep 
imaging observations from ~2 pm to millimeter wavelengths, 
including the submillimeter domain, are required. 

The Herschel Space Observatory now makes it possible to 
bridge the submillimeter gap by providing deep, wide-field dust 
continuum imaging at 6 wavelengths from 70pm to 500pm. In 
particular, the Herschel Gould Belt survey ( Andre et al.|2010| is 
imaging the bulk of the nearest (d < 0.5 kpc) molecular cloud 
complexes with the SPIRE and PACS cameras and will eventu- 
ally provide a complete census of Class protostars and prestel- 
lar cores in these clouds. 

Here, we present a millimeter continuum study of the Aquila 
Rift complex which addresses the question of the lifetime of the 
Class phase and that of the star formation rate in two nearby 
cluster-forming clumps: the Serpens South and W40 protoclus- 
ters. 

The Aquila Rift complex 

The Aquila Rift molecular cloud complex corresponds to a large 
extinction feature (see [Bontemp s et al.|[2010| l, located above 
the Galactic plane at 1=28°, only ~3° south of the well-known 
Serpens Main star-forming region. The Aquila Rift extinction- 
defined area (see Bontemps et al. 2Q10| > mainly harbours two 
known sites of star formation: Serpens South is the western 
young embedded cluster revealed by infrared Spitzer observa- 
tions (Guterrr mfh et al.||2008[ ), and W40 is the eastern clus- 
ter associated to the eponymous Hn region, also known as 
Sharpless 2-64 ( [Smith et al . 1 1 9851 1 Vallee 1 1 987] > . 

A recent Herschel mapping of ~ a 3.3°x3.3° field was carried 



out toward this complex as part of the Gould Belt survey (Andre 
|et al.|201Q] >. The corresponding far-infrared maps have unprece- 
dented spatial resolution and sensitivity in 5 bands between 70 
and 500 pm. The first results of these Herschel observations of 
the Aquila complex are presented in Bontemp s et al.| ( |2010| l, 



and Konyves et al. ( 2010] >. Briefly, the Herschel maps revealed 



a large number of new protostars and prestellar cores mainly 
spread over three currently active star-forming sites: the filamen- 
tary molecular cloud harboring the Serpens South infrared em- 
bedded cluster in its center, the eastern HII region W40/Sh2-64, 
and the southernmost MWC297 / Sh2-62 region. 



Distance of the star forming clouds in the Aquila Rift 

The distances of the Serpens South and W40 protoclusters, both 

located inside the extinction feature of the Aquila Rift, are still 

uncertain. 

Whether or not the Aquila Rift complex and the Serpens 
Main cloud are part of the same complex and therefore at the 
same distance is still a matter of debate ( |Gutermuth et al.|2 008 ; 
Dzib et al.|2010|[Loinard et al.|201 1 1, however from the extinc- 



tion maps obtained by Bontemps et al. ( 2010] >, their respective 



extinction features are seen as clearly distinguishable separate 
regions. Indeed, assuming that Serpens Main and Aquila Rift 
are different entities also has the advantage of explaining the 
discrepancy found between their extinction-based distance es- 
timates ( |Straizys et al.|1996| [20031 |Bontemps et al.|2010"l|Knude 
[2011) : -200-280 pc for the Aquila Rift complex (due to the 
known extinction wall at a distance of 225 ± 55 pc), and ~415 pc 
for Serpens Main (from the recent VLBI-based distance by Dzib 
|et al.|2010 1, suggesting that Serpens Main is located behind the 
extinction wall associated to Serpens South. 

Previous studies targeting the Serpens South protocluster 
(Gutermuth et al. 2008 Bontemps et al. 2010 Konyves et al. 
2010) used the extinction-related distance of ~260 pc for this 
star forming region, Serpens South being the highest extinction 
region found inside the Aquila Rift complex. Also based on ex- 



tinction measure ments (|Smith et al. 1985[ l and radio spectral line 
data ( |Vallee|1987[ | Vallee et al.|1992) , the nearby HII region W40 
has usually been considered to be at a distance ranging from 300 
to 700 pc. A recent Chandra study by |Kuhn et al. ( 2010] > em- 
phasizes, on the basis of X-ray luminosity function arguments, 
a most-likely distance of ~600 pc for the W40 cluster, although 
not excluding the possibility of a smaller (~300 pc) distance. 

We argue that both Serpens South and W40 are likely to be 
at the same distance because (i) they belong to the same continu- 
ous extinction feature (see Bontemps et al.|20T0"l|Schneider et alT 
201 1 ) associated to the Aquila Rift, and (ii) they are populating 
a continuum of LSR velocities between 4 and 10 km.s -1 . 



For consistency with the Herschel Gould Belt survey (Andre 



et al. 2010 1 results for that region (Bontemps et al. 2010 
Konyves et al.|2010 1 and easier comparison with previous stud 



ies, we adopt a distance of 260 pc for the Aquila Rift complex, 
including both Serpens South and W40, throughout this paper. 



2. MAMBO 1.2 mm dust continuum observations 
and data reduction 

We used the Max-Planck Millimeter Bolometer array 
(MAMBO) on the IRAM 30-m telescope near Pico Veleta, 
Spain) during March 2007 to carry out a millimeter continuum 
mapping of the Aquila Rift complex. Ten adjacent sub-fields 
were selected thanks to the near-IR extinction map derived 
from 2MASS data ( |Bontemps et al.|[20T0l |Schneider et al. 
201 1| ), in order to image the dust continuum associated with 
the main extinction features of the Aquila Rift complex, at 
1.2 mm (~240 GHz). These sub-fields were mapped using the 
on-the-fiy technique in the double-beam switching mode, i.e. 
chopping was performed by switching the wobbling secondary 
with a throw ranging from 62" to 160". The scanning speed was 
8"/sec. The atmospheric optical depth at zenith was estimated 
between 0.09 and 0.41 at A = 1.2 mm. Pointing and focus 
checks were made every 3 hours, resulting in a pointing accu- 
racy better than 3", and the calibration obtained by mapping 
Mars (primary calibrator) results in a flux accuracy better than 
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Right ascension (J2000) 

Fig. 1. 1.2 mm dust continuum MAMBO mosaic of the Aquila Rift complex, smoothed to an effective angular resolution of 13". 
Due to the high dynamic range of fluxes in the map, and to enhance visibility of the structure of emission in the whole area 
mapped, the color scale chosen here has a maximum of 400 mJy/beam. The highest flux density of order ~1400 mJy/beam is found 
in the center of the Serpens South protocluster (see the blow-up in Figure [2] In this smoothed map, the mean rms noise level is 
~ 13 mJy/13"-beam. 



~ 15%. The angular resolution with IRAM 30-m antenna is ~ 
11" (HPBW) at A = 1.2 mm, and the average rms achieved is 
~15 mJy/1 l"-beam. We reduced and co-added all the MAMBO 
sub-fields using standard methods in MOPSIC, part of the 
GILDASQ software, and produced a global mosaic covering a 
total area of ~ 60' x 50' (see Fig.[T]). 



3. Dust continuum mapping: results and analysis 

3.1. Main features of the millimeter continuum maps 

The large-scale 1.2 mm mosaic obtained with the MAMBO cam- 
era is shown in Fig.[T] Four main parsec-scale regions of interest 
can be identified in this mosaic, which are labelled 1,2, 3, and 4. 
Indeed, these regions show both extended emission and compact 
structures in dust continuum, and are therefore likely tracing ac- 
tive star-forming areas. 

Region labelled 1 is the previously known Serpens South 
protocluster, showing up as an elongated ~1 pc filament in the 
1.2 mm dust continuum map (see blow-up in Fig. [2]), coinci- 
dent with the Spitzer near-infrared absorption feature discussed 
by Gutermuth et al.| ( |2008[ >. Very condensed in its center, the 
filament sees its density decreasing while going along the two 
arms that have a cross-shape in more diffuse emission. Compact 
sources are found up to 15' away from the center of the filament. 



at the convergence zone of the filament shows several dust con- 
tinuum peaks at the MAMBO resolution: it is the center of a pro- 
tocluster, with several compact millimeter sources tracing very 
young protostars at the beginning of their accretion phase (sepa- 
rated by less than 5000 AU), as well as numerous more evolved 
infrared YSOs. 

Region labelled 2 is shown with much greater details in 
Fig. [3] It mainly features a ring-like structure following the outer 



border of the W4 HII region ( |Smith et al.|[1985] |Valleep987 
Kuhn et aT7||2010|>, and spatially coincident with the CO emis 



sion mapped by |Zhu et al.| ( |2006[ l. In addition to this western 
dense shell, the map shows elongated structures (also contain- 
ing compact sources) in the eastern part of the photo-dissociated 
region, while the central part of the infrared cluster containing 
most of the pre-main sequence objects ( |Kuhn et al.|20 10) does 
not show any evidence for cold dust emission at 1 .2 mm. 

Finally, regions labelled 3 and 4 are more diffuse, less rich 
star-formation sites, for which no previous literature was found. 

The ~ 1 l"angular resolution of MAMBO at 1 .2 mm allowed 
us to resolve individual condensations within the dust cores at 
the assumed distance of Aquila Rift, and we therefore used our 
millimeter data to complement the Herschel dataset and charac- 
terize the protostellar population in the mapped area. 



This high-density region (column densities up to 3xl0 23 cnr 2 ) 3.2. Source extraction and properties 



1 Grenoble Image and Line Data Analysis System, software provided 
and actively developed by IRAM (http://www.iram.fr/IRAMFR/ 
|GILDAS| l 



In order to improve the detection of compact sources in our 
field, we smoothed the original 1.2 mm continuum map to 
a 13"-beam, and performed a systematic source detection on 
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Right ascension (J2000) 

Fig. 2. Blow-up of the Region 1 from the MAMBO 1 .2 mm dust 
continuum map shown in Fig.[T| around the Serpens South fila- 
mentary cloud. The mean rms noise level is ~ 12 mJy/13"-beam. 
Black contours show the levels of 5<x, 8<r, 15cr, 30cr, 70cr, 90cr, 
and 120cr (1.44 Jy/beam). 




18 31 30:0 
Right ascension (J2000) 

Fig. 3. Blow-up of the Region 2 from the 1.2 mm continuum 
map shown in Fig.[T] around the HII region W40. The mean rms 
noise level is ~ 14 mJy/13"-beam. Black contours show the lev- 
els of 5cr, 8cr, 15cr, 20cr, 25<r, and 35<x. One of the candidate 
exciting stars, Sh64-OS2a, is detected in our MAMBO mapping 
and indicated here. 



the MAMBO map using the new source e xtraction method 
getsources PIl Men'shchikov et al. |26To| |20 1 1 ) > . This analysis led 



2 This new algorithm analyzes fine spatial decompositions of origi- 
nal images across a wide range of scales and across all wavelengths, 
cleans those decomposed images of noise and background fluctuations 
separately at each wavelength, and constructs a set of wavelength- 
independent detection images that preserve information in both spa- 
tial and wavelength dimensions as well as possible. Sources are de- 
tected in the set of decomposed combined detection images by follow- 
ing the evolution of their segmentation masks across all spatial scales. 
Measurements of the sources' properties are performed in the original 
images at each wavelength after the background has been subtracted 



to the extraction of a total of 77 dust continuum sources in the 
whole region shown in Fig.[T] among which 25 are located in the 
Serpens South filament, while 36 lie in the W40 region. The re- 
sulting MAMBO compact sources catalog is reported in Table 1 . 
Starting from this 1.2 mm dust continuum set of sources, par- 
allel measurements were performed at other wavelengths with 
getsources in order to estimate the counterparts of the MAMBO 
sources in both the Herschel (from 70^m to 500yum), the Spitzer 
(8yum and 24//m) and the 2MASS K-band maps, when data was 
available. All these sources were also examined individually by 
eye, ensuring association to the right counterparts, which is not 
always straightforward while working with different sizes of 
beams in crowded regions, such as Serpens South or W40. 

We used these flux measurements to establish spectral en- 
ergy distributions (SEDs) for all the 77 MAMBO sources, and 
derived other source properties (dust temperature, core/envelope 
mass, bolometric temperature and luminosity) when a good fit 
was found to reproduce the SED obtained between 2/itn and 
1.2 mm with grey-body functions. The condensation masses 
(core masses for starless sources and envelope masses - M env 
- for protostellar objects) and dust temperatures Td us t were es- 
timated from the fitting of these SEDs by grey-body func- 
tions (Ko nyves et al.|20T0]>. W e adopted the recommended val- 
ues (Ossenkopf & Henning 1994 Henning et al. \995) for 
the dust opacity per unit mass column density (of dust and 
gas): Ki.2mm = 0.01 cm 2 g _1 for protostellar envelopes, and 
Ki.2mm = 0.005 cm 2 g _1 for starless clumps. The condensations 
have masses ranging from ~ 0.02 M Q to ~ 6 M , with an aver- 
age factor ~ 2 absolute uncertainty (see Table 1). Note also that 
for the most evolved sources (mainly Class I protostars), these 
masses might be contaminated or even dominated by disk mate- 
rial, while they represent the core masses for the starless sources. 

The MAMBO sources showing no counterpart shortward 
of 160 yum with Herschel-PACS, but detected in at least two 
Herschel- SPIRE bands (250, 350, or 500//m) were classified 
as starless sources. The MAMBO sources with compact coun- 
terparts at Herschel wavelengths < 160 fj.m were classified as 
candidate young stellar objects (YSOs), The bolometric lumi- 
nosities Lboi, submillimeter luminosities L^^ o 0fl , and bolometric 
temperatures T\, \ of the candidate YSOs were estimated by in- 
tegrating their SEDs, and are reported in Table 1 . Note that the 



values of L 



and Lboi derived here are expected to be more 



accurate (estimated uncertainty <40%) than previous estimates 
from the literature, thanks to the Herschel data providing a good 
coverage of the far-infrared/submillimeter domain. 

3.3. Evolutionary classification of protostellar sources 

All of the candidate embedded protostellar objects detected with 
MAMBO were then placed in the envelope mass versus bolo- 
metric luminosity (M env -Lboi) diagram shown in Fig.|4] in or- 
der to distinguish between envelope-dominated (M env >> M*) 
Class protostars and more evolved (Class I) protostars. In this 
diagram, the bolometric luminosity of a protostellar object can 
be used as an indirect tracer of the stellar mass M* it has al- 
ready accreted ( ] Andre et al|l993| |Andre & Montmerle||1994| l, 
based on comparison with simple accretion model tracks (cf. 



Bontemp s~et al.|1996| |Saraceno et al.|1996| | Andre et al.pO OO). 
The specific evolutionary tracks shown in Fig.|4|are similar to the 
tracks used by Andre et al. ( 2008 ) and were computed assuming 
each protostar forms from a prestellar core of finite initial mass 



by interpolation under the sources "footprints" and overlapping sources 
have been deblended, see Men'shchikov et al. 201 1 (in prep). 
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Fig. 4. Envelope mass versus bolometric luminosity diagram comparing the locations of the protostars detected with MAMBO in 
the Aquila Rift region (open triangles for the Serpens South sub-field, open squares for the other sub-fields) with the positions of 
knwon low-mass Class I (filled circles) and Class objects (open circles) in Ophiuchus, Perseus, and Orion ( [Andre et al.|2000) . 
Numbers reported inside triangles are the reference numbers of the sources reported in the last column of Table 1 . Model protostellar 
evolutionary tracks, computed for various final stellar masses (indicated above the Lboi axis) assuming the accretion/luminosity 
history described in the text, are superimposed (4 blue tracks representing protostellar evolution for objects with final masses 0.1, 
0.3, 1 and 3M ). Evolution proceeds from the upper left to the lower right as indicated by two arrows on each track, plotted when 
50% and 90% of the final stellar mass has been accreted, respectively. The straight lines show two M em — Lb i relations marking the 
conceptual border zone between the Class (M env > M*/e) and the Class I (M env < M*/e) stage, where e is the local star formation 
efficiency taken here to be 50%. The dashed line is such that M em , oc L bol (cf. Andre et al.|1993| Andre & Montmerle 1994 ]), whil e 
the dotted relation follows M env oc as suggested by the accretion scenario adopted in the tracks (see also Bontemps et al.|l996| . 
The typical uncertainty in the plotted values of Lboi is ~ 30 - 40%, while the uncertainty in the M env values is a tactor ot ~ 2 (see 
§ 3.2). The black arrow in the upper right indicates how the positions of the W40 protostars would move in the diagram if a distance 
of 600 pc was adopted for W40, instead of 260 pc. 



M core = M env (0) and has L bo , = GM*(f)M acc (f)/fi*(f) + L+(f), 
where /?+(M*) is the protostellar radius and L* the PMS birth- 
line luminosity (Stahler T988J. The mass accretion rate and 
the envelope mass were assumed to be related by M acc (f) = 
£ M em (t)/r, where e = 50% is the typical star formation effi- 
ciency for individual cores (cf. Matzner & McKee 2000) and 
t - 10 5 yr is the characteristic timescale of protostellar evolu- 
tion, leading to M acc (f) and M env (f) functions declining exponen- 
tially with time as suggested by the observed decline of outflow 
power from Class to Class I (see Bontem ps et al.|19 96 ). Note 
that we adopt e = 50% here a s in|Andre et al.| ( |2008| ) instead of 



sified as Class I or Class II, depending on their infrared spec- 
tral index am = dlogAF^/dlogA between 2.2/im and 24/^m (as 
estimated by us from independent getsources measurements at 
2.2/im, 8/im, 24/mi when available - see §4). 
The use of this evolutionary diagram (see Fig. [4]) leads to the 
identification of 9 Class (and 3 border-line Class 0/1) proto- 
stars in the Serpens South filament. Likewise, 8 Class (and 4 
border-line Class 0/1) protostars are found in the W40 region. 
The other sub-fields labelled 3 and 4 gather 8 Class (and 2 
border-line Class 0/1) protostars. 



e = 100% as in [Andre et al.| ( |2000| ), since recent comparisons |2000b and [7j 



Two other criteria, [L^^/L ho i 



of the prestellar core mass function with the stellar IMF suggest 
e<50% (e.g. l Alves et al.|2007j | Nutter & Ward-Thompson 12007 1 
Konyves et aL|2UlU||Andre et al.|2010| ). 



bol 



> 1 

70 K] "(e.g. [Myers et 



%f1 
al. T 



I (e.g. 



Andre et al. 



998) respectively, 



A Class status was attributed to sources with M env > M*/e 
according to the diagram of Fig.|4] i.e., protostars which have 
accreted less than half of their final stellar mass. Conceptually, 
protostars at the border line between the Class and the Class I 
stage are objects which have accreted exactly half of their final 
stellar mass, and which thus have M env = M*/e according to 
the tracks. In practice, sources falling into the "Border Zone" 
region located between the dashed and dotted lines in Fig. [4] 
have been classified as border-line Class 0/1 objects in Table 1. 
Sources lying below this border zone in the diagram were clas- 



have been commonly used in the literature to select candidate 
Class protostars among embedded YSOs. In order to assess the 
importance of uncertainties in our classification, we have applied 
these two alternative methods of classifying embedded YSOs to 
the Aquila Rift sample and compared the results to those ob- 
tained above on the basis of the M env - Lboi diagram. 
We find (see Table 1) that using the L^^/L bo i > 1% crite- 



The original criterion introduced by 



L ^bmm / L boi > 0.5%, but only YSOs with L^ h "^/L bo i > 1% were con- 
sidered bona-fide Class objects by Andre et al. (2000), while sources 
with L^7L b0] between ~ 0.5% and ~ 1% were classified as border- 
line objects (see Table 1 of Andr e et al.|2000| . 



Andre et al 
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Fig. 5. Same as Figure 4, but showing only the sources from the Serpens South cluster forming clump. Triangles show the locations 
of the MAMBO protostars (green triangles are Class and red triangles are Class I protostars, while black triangles represent the 
border-line objects). Red circles show the locations of the Class I protostars detected by Spitzer and without MAMBO detection. 
The Class I sources without any MAMBO counterpart are indicated as upper-limit values for M env , around the detection limit of 
2.5xlO~ 2 M . The green double arrow along the M env axis shows the range of masses for the prestellar cores detected by Herschel, 
based on a preliminary analysis of the Gould Belt survey data, which are 85% complete down to 0.3M o in this region (Konyves 
et al.|2010) . The respective locations of the prestellar cores, Class and Class I protostars in the Serpens South protocluster are 
roughly consistent with what is expected from a single population of objects evolving along the model evolutionary tracks shown in 
blue. 



rion (which traces a ratio of envelope to stellar mass larger than 
unity - see Andre et al.| 1993 1 to select Class protostars leads 
to similar results to those obtained from the M env - Lboi diagram. 
Only in a few cases does the use of this criterion lead to the 
selection of border-line Class 0/1 sources (according to M env - 
Lboi diagram) as candidate Class protostars. For example, us- 
ing the L > ^ 5 ° M / L^i > 1% criterion leads to the selection of 12 
Class protostars in the Serpens South filament (compared to 9 
Class and 3 border-line sources when using the M env - Lbd di- 
agram), and 1 1 Class protostars in the W40 region (compared 
to 8 Class and 4 border-line sources when using the M env - Lboi 
diagram). Based on the present study, we suggest that setting the 
border line at L^^/Lboi > 3% may better match the results 
obtained with the more conceptual classification based on the 
tracks in the M em - Lboi diagram that we use here. 
We find that using the Lboi < 70 K criterion for distinguishing 
Class protostars from Class I objects ( |Chen et al.| 19 95 ; Myers 
|et al.|1998] ) tends to systematically overestimate the number of 
candidate Class objects, compared to the numbers obtained 

from both the M env - Lboi diagram and the / Lboi > 1 % cri- 

terion. This is especially true in the "warm" W40 region where 
the estimates of Lboi are more uncertain due to confusion at mid- 
infrared wavelengths from the photo-dissociation region. For ex- 
ample, using the Lbd < 70 K criterion leads to the identification 
of 12 Class protostars in the Serpens South filament, and 19 
Class protostars in the W40 protocluster (twice the number of 
Class protostars found in this region if using the two other cri- 
teria presented above). 

We therefore conclude, in agreement with |Young & Evans| 
(2005), that the observational criterion to distinguish Class 



protostars making use of the L x> ? 11 /Lboi ratio seems to be more 
reliable than the 7boi criterion, also because it is less sensitive 
to environment and geometry effects. Comparing the various 
methods of classifying protostars, which is roughly equivalent 
to taking into account the border-line Class 0/1 protostars from 
the M env - Lbd diagram, we find that the typical uncertainty in 
the number of Class protostars is ~50% in W40, and ~30% in 
Serpens South. 



4. Discussion and conclusions 

4.1. Class phase lifetime and star formation rates 

Here we propose a statistical analysis of the properties of the 
protostellar populations in the two main regions showing active 
star formation from our MAMBO map: the Serpens South fila- 
mentary cloud and the W40 HII region. 

4.1.1. The Serpens South protocluster 

First, we take advantage of having a complete population of 
Class protostars in a spatially well-defined star-forming cloud 
to derive an estimate of the lifetime of the Class phase. 

The Serpens South protocluster was observed with IRAC 
and MIPS onboard the Spitzer telescope, and the analysis of 
IRAC data by |Gutermuth et al.| ( |2008] l revealed 37 Class II YSOs 
and 54 Class I protostars (including flat spectrum sources) in 
a region of 14'xl0' around the filament. Based on the pub- 
licly available MIPS-24jt/m, IRAC-8jum, and 2MASS//T-band 
data, we independently performed an extraction of MIPS-24//m 
point-like sources with getsources in the 0.8 x 1.2 pc area 
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around the Serpens South filament shown in Fig.|2j which re- 
vealed a total of 112 objects with an infrared spectral index 
am = dlogAFj/dlogA > 0.3 between 2/im and 24pm, that we 
considered as candidate Class I protostars. 
As discussed in § 3.2 above, our MAMBO mapping finds 9 can- 
didate Class protostars in the same area (not counting border- 
line Class 0/Class I sources). Our census of Class protostars in 
Serpens South is sensitive to protostars with final masses down 
to ~ 0.1 M G (see evolutionary t racks in Fig. [4| . According to 
the grid of protostellar models by Robitaille et al. (2007 ), Class I 
protostars with final masses M^ nal > 0.1 M Q are expected to have 
24/itn flux densities F^ifoa > 35 mJy at d — 260 pc (assuming 
a random distribution of inclination angles, i.e computed using 
a mean value of ~57°). We will therefore consider here only the 
48±2 candidate Class I sources which have F24f, m > 30-40 mJy 
from our extracted sample. Note that at such a 24//m flux den- 
sity level, our Class I sample should not be contaminated by 
extragalactic sources (cf. Harvey et al.|2007 i. Among these 48 
objects, 16 coincide with 21 IRAC Class I sources reported by 
Gutermuth et al. (2008 ), while 21 are new detections made pos- 
sible thanks to the use of the 24 pm data. The Spitzer, Herschel, 
and MAMBO maps are very likely wide enough to probe com- 
plete samples of YSOs associated with the Serpens South fila- 
mentary clump, since there are no YSOs, neither high column- 
density gas, immediately outside this 0.8 x 1 .2 pc area (shown in 
Fig. [2} in available larger Herschel, Spitzer, and MAMBO maps 
of the region. 

Assuming that the star formation rate in Serpens-South has been 
constant over the duration of the protostellar phase, we can esti- 
mate the lifetime of the Class phase by comparing the relative 
numbers of Class and Class I sources in this region. Combining 
the MAMBO and Spitzer samples discussed here, the Serpens 
South cluster-forming filament harbors a total of 57 protostellar 
objects with expected final stellar masses > 0.1 M , and more 
precisely 84% of Class I and 16% of Class objects. Adopting 



a lifetime of 0.2-0.5 Myr for the Class I phase ( |Greene et al. 
1994"} |Evans et al.]|2009"] ), the ratio of Class to Class I proto- 
stars we find, N(Q)/N(I) = 0.19, suggests a typical lifetime of 
4-9 xlO 4 yr for the Class phase. This r esult is intermediate 



between the early estimate of 1-3 xlO 4 yr by |Andre et al. (2000) 



and the more recent estimate of 0.10-0.16 Myr by [Evans et al. 
( |2009l ) 

Our analysis makes use of an unprecedented census of YSOs 
in Serpens South, using wavelengths from 2 pm to 1.2 mm and 
including the Herschel submillimeter domain. It should there- 
fore provide a more complete population of embedded YSOs in 
Serpens-South and a more reliable assessment of their evolution- 
ary states (via better determination of their envelope masses and 
bolometric luminosities) than earlier studies. 
However, the lifetime we derive for the Class phase is still 
subject to large uncertainties, partly due to remaining uncertain- 
ties on the lifetime of the Class I phase. Estimates of the latter 
are themselves highly dependent on assumptions concerning the 
lifetime of the Class II (or T Tauri) phase (see |Evans et al. 2009 ). 
The method chosen to distinguish Class protostars from more 
evolved YSOs also impacts the determination of this lifetime, 
because different methods lead to slightly different source counts 
(see § 3.3). However, the two favored methods to classify YSOs 
were compared in § 3.3, and the differences they show (typical 
~30% uncertainty on the number of Class protostars detected, 
i.e N(0)/N(l) ~ 0.19 - 0.27) are expected to result in a typical 
uncertainty of ~0.02 Myr in the Class lifetime we find here, 
which is a factor of 2 lower than the main source of error due to 
the uncertainty in the Class I lifetime (leading to a ~0.05 Myr 



uncertainty on the Class lifetime as shown above). Moreover, 
we do not claim our result to be universal, as our study is fo- 
cused on only one star-forming region and does not address the 
question of environmental effects. 

Even taking these caveats into account, we stress that our results 
tend to support a scenario of relatively fast accretion at the be- 
ginning of the protostellar phase in protoclusters (cf. |Henriksen| 
|et al.||1997] ), since they suggest that a typical YSO accretes the 
first half of its final stellar mass in < 0. 1 Myr (Class phase) 
and the second half in ~ 0.4 Myr (Class I phase - cf. Evans et al. 



|2009| >. 

Finally, we remark, in Fig. [5] an over-abundance of low- 
luminosity Class I protostars if compared to what is expected 
from the population of Class protostars and associated model 
evolutionary tracks in the M env -Lboi diagram for Serpens South, 
which is likely to be an illustration of the well-known "lumi- 
nosity problem" first noted by |Kenyon et al.| ( fl990[ ) and more re- 
cently stressed by Evans et al. ( 2009 ). However, we point out that 
the locations of prestellar cores, Class and Class I protostars in 
the M env -iboi diagram for Serpens South shown in Fig. 5, are 
roughly consistent with what expected for a single population of 
objects evolving along the model evolutionary tracks. This illus- 
trates the validity of our method for determining the evolutionary 
status of protostars. 



Second, using the column-density map of the Aquila Rift 
derived from Herschel data (Konyves et al.pOlO i 



we can es- 
timate a mass for the Serpens South filamentary cluster-forming 
clump shown in Fig. [2] and derive a value of the star for- 
mation efficiency (SFE) for this region. We find a total mass 
of ~610 M Q in the 15' x 25' area containing the protostel- 
lar population discussed here, with a mean column density of 
Nh 2 ~ 3. 1 x 10 22 cirT 2 . Considering that this region has formed 9 
Class 0, 48 Class I, and 37 Class II objects, i.e. a total of 94 YSOs 
of mean mass 0.5 M , we find: S¥E-M slar l{M star + M,, as ) = 
47 / (47 + 6 1 0) ~ 7% . Assuming a standard duration of 2 Myr for 
going from the Class to the Class II phase ( |Evans et al.|2 009 ) 
and a projected area of ~ 1 pc 2 for the region harboring YSOs, we 
derive the following star formation rate per unit surface for the 
Serpens South cluster-forming clump: SFR~23 M Q Myr~'pc" 2 . 
This value is significantly higher than the typical SFR values 
reported by Evans et al. (2009) for entire molecular clouds 
(SFR~1.6 M Q Myr^pc"" 1 ), which is consistent with the fact that 
we are focusing here on a compact cluster-forming clump with 
higher column densities (up to Nh^ ~ 10 23 cm 2 ) than the clouds 
targeted by |Evans et al.| ( |2009[ ), also showing evidence of a re- 
cent burst of star formation. A recent burst of star formation is 
indicated by the unusually large number of Class and Class I 
protostars found in the Serpens South star-forming clump com- 
pared to the relatively small number of Class II YSOs, which 
was already pointed out by |Gutermuth et al. (2008) and that we 
confirm here. Note that if Serpens-South has indeed undergone a 
recent burst of star formation, then the observed ratio of Class 
to Class I protostars (N(Q)/N(I)) may overestimate the actual ra- 
tio of Class to Class I lifetimes. 

Finally, the SFE value derived here for Serpens South is of the 
same order as that reported by Evans et al. ( 2009[ ) for molecu- 
lar clouds (i.e. SFE~3-6%). The fact that Serpens South shows 
both a high SFR and a relatively low SFE may be interpreted as 



an additional sign of its very early evolutionary stage (Lada & 
|Lada|2003| ). 
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4.1.2. The W40 Hll region 

The present MAMBO 1.2 mm continuum study allows us to 
identify 8 candidate Class protostars in the W40 cluster region, 
where copious photo-dissociation processes (Hll region proba- 
bly powered by the massive young star Sh64-OS2a, see Pig. [3]» 
make it difficult to obtain a census of YSOs at shorter wave- 
lengths. Indeed, the Lyman-continuum luminosity towards the 
Hll region of W40 (~1.5xl0 48 photons.s" 1 ) inferred from the 
radio continuum emission by Altenho ff et al.|(fl 970) could be 
produced by a single 09 V star (Smit het al.|1985[ Tln a recent 
study, |Kuhn etal] ((20T0|> found 93 YSOs with NIR ZT-band ex- 



cess belonging to the W40 region, based on their Chandra ob- 
servations toward the young cluster. However, their census did 
not distinguish between Class I, Class II, and later evolutionary 
stages. It is therefore difficult to derive accurate lifetimes for the 
Class and/or Class I phases in the W40 region. 

We estimate a gas mass of 310 M Q from the Herschel col- 
umn density map ( |Konyves et al.|[20T0| of the region shown in 
Fig.|] s imilar to what wa s found (320 M ) from the 13 CO(2-l) 
map by |Zhu et al. ( 2006| l. The mean column density over this 
region is Nh 2 ~ 1-8 X 10 22 cm~ 2 . Considering that this clus- 
ter presently contains a total of 20 protostellar sources detected 
with MAMBO and 93 Class II so urces showing NIR K-baxvd 
excess and detected with Chandra fKuhn et al.|2010| l, we esti- 
mate that the current star formation efficiency of the protocluster 
is SFE = 56/(310+56) ~ 15% assuming a mean stellar mass of 
0.5 M Q ^\ We stress here that the mass derived from the Herschel 
submillimeter continuum maps may be significantly lower than 
the original clump mass reservoir that formed the PMS objects 
presently observed, as a fairly large fraction of the gas may have 
already been expelled from the region. This is suggested for ex- 
ample by the lack of any submillimeter/millimeter continuum 
emission toward the center of the cluster where most of the PMS 
stars are detected. 

Assuming a distance of 260 pc for W40, the projected area 
of the W40 region harboring YSOs (area shown in Fig. [3j is 
0.9 pc 2 . Assuming a period for star formation of 2 Myr (up to 
the Class II population), we can derive a lower-limit value for 
the SFR of 34 M Myr _1 pc~ 2 , which is similar to what is found 



for the Serpens South protocluster (see §4.1.1). However, Kuhn 
et al. (2010 1 argue that the high disk fraction they find thanks to 



/T-band excess measurements would support an age for the W40 
cluster of only ~ 1 Myr, which would lead to a upper-limit SFR 4 
of -57 M Q Myr -1 pc~ 2 . For comparison, the star for mation rate in 
the Serpens Main cluster, i.e. 56 M Myr 1 pc -2 (Harvey et al. 
2007 ), is of the same order of magnitude as the values we find 
here. 

We note that the classification we adopted here to determine 
the evolutionary status of the MAMBO sources, based on the 
M env - Lboi diagram shown in Fig. [4] is independent of the dis- 
tance used for W40 as the envelope mass and bolometric lumi- 
nosity have the same dependency with distance. At first sight, 
we may expect the SFR and SFE values estimated for the W40 
complex to scale as I Id 2 , since both the gas mass and the cloud 
area scale as d 2 , while the number of YSOs detected stays the 
same. However, assuming a standard IMF, the total stellar mass 
is roughly proportional to the total bolometric luminosity of the 
YSO population, which also scales as d 2 . Therefore, we argue 



4 Note that if we take into account the additional Class III objects 
identified in X-rays by |Kuhn et aLlpOTO) (-109 Chandra YSOs with- 
out A^-band excess ), then the SFE goes up to -26% and the upper- 
limit of the SFR to ~120 Mq Myr~'pc~ 2 (assuming a lower-limit age of 
~1 Myr, see [Kuhn et al.|20To| . 



that both the SFR and the SFE we estimate for W40 do not de- 
pend much on the adopted distance. 



4.2. Comparing evolutionary models and IMF predictions in 
our Serpens South sample 

Considering that we detected 57 protostellar objects with final 
masses > 0.1 M in the Serpens South cluster-forming 
clump, and following the IMF ( |Kroupa|[200T] |Chabrier||2003] l, 
one would expect the most massive object from our sample to 
have M^™^. ~ 2.5±0.5 M (taking into account the -30% un- 
certainty on the number of protostellar objects in Serpens South). 
Similarly, applying the IMF to our sample of YSOs predicts a 

median final mass M^ nal - 0.35±0.1 M , for our sample. 
The masses of most of the Class I sources are not constrained by 
our MAMBO observations (43 out of 48 are Spitzer sources not 
detected with MAMBO). Moreover, variable accretion rates dur- 
ing the embedded phases (only scenario likely to solve the lumi- 
nosity problem, see §3.3, seen in our data by an over-abundance 
of low-luminosity Class I protostars in Fig. [5]), do not allow us to 
use luminosities to derive accurate masses for Class I protostars 
not being detected at sub-millimeter or millimeter wavelengths. 
Therefore, we use the properties from our sample of Class pro- 
tostars to compute the median stellar mass and maximum stellar 
mass, expected from the model evolutionary tracks of Fig. [5] We 
then compare these values with values expected from the stan- 
dard IMF. We stress however that the low number of Class 
protostars we use here is not allowing to draw robust conclu- 
sions. 

We find a median envelope mass M env ~ 0.8 M . If we follow 
the evolutionary model tracks shown in Fig. [5] we predict a final 

median mass of M^ nal ~ 0.4 M . Moreover, the most massive 
Class object from our sample (SerpS-MM18) has M™ x -3+1 
M , and evolutionary models therefore suggest a maximum final 
stellar mass -1.7+0.5 M . 

Thus, the use of the model evolutionary tracks drawn in Fig. [5] 
predicts both a median and maximum stellar mass in agreement 
with a random sampling from the standard IMF ( |Kroupa|2 001 ; 
|Chabrier|2003l . 

4.3. Two modes of efficient clustered star formation? 

The high values of SFRs we derive for both Serpens South and 
W40 are consistent with the fact that we consider dense and com- 
pact regions actively forming large numbers of protostars from 
large mass reservoirs, i.e protoclusters (cf. Hurt & Barsony 1996; 



Mye rs|2001| |Andre et al.|2000} |Lada & Lada|2003| l. Theselwo 
adjacent star-forming regions are likely to share the same dis- 
tance and have similar projected areas with ongoing star forma- 
tion as well as similar gas reservoirs (gas masses of 610 M and 
310 M , respectively). However, the formation of each of these 
two young protostellar populations is likely to have been initi- 
ated in a different way: 

(1) Large-scale collapse of a gravitationally unstable filament 
for the Serpens South cluster-forming clump, similar to what 
was observed in the intermediate-mass protocluster NGC 2264- 
C ( |Peretto et al.||2006| >. Indeed, it has been proposed based on 
the first results from Herschel ( |Andre et al.||2010| ) that a dense 
filament with a highly supercritical mass per unit length (cf. 
Inutsuka & Miyama 1997) such as Serpens South may become 
globally gravitationally unstable and undergo overall collapse, 
eventually leading to the formation of a rich protocluster. Our 
MAMBO 1.2 mm map confirm that the high-column-density fil- 
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amentary clump surrounding the central mid-infrared embedded 
cluster detected by Spitzer ( Gutermuth et al.|2008 1 in Serpens- 
South is actively forming young Class protostars. The high 
fraction of Class and Class I protostars compared to the to- 
tal number (30-40) of Class II YSOs, found in the central core 
of the protocluster but also in its filamentary dusty arms, is 
suggestive of a recent burst of star formation in this cluster- 
forming clump. This view is also supported by the recent results 
of Nakamura et al. (201 lb ), who detected numerous protostellar 
outflows toward this clump. 

(2) Pressure-triggered star cluster formation (see [Elmegreen 



1998] |Deharveng eTaL|2005] > around the HII region W40. The 



winds and energetic radiation from the massive young stars of 
W40 have created an expanding HII region and compressed 
the gas clump at the origin of the cluster, forming a ring-like 
structure of dense gas now harboring a second generation of 
young protostars. Indeed, all of the youngest (Class 0) YSOs 
are found in the high column-density regions surrounding the 
infrared cluster containing most of the more evolved Class II 
and Class III objects, while not many of these more evolved ob- 
jects are found in the outer dense regions harboring the Class 
protostars. If, as suggested by the high disk fraction reported by 
Kuhn et al. (2010 1, the young infrared cluster is indeed ~1 Myr 



old, then W40 is currently undergoing a small burst of star for- 
mation, and may end up forming a rich bound cluster. 

We have shown that these two protoclusters, although re- 
sulting from different collapse initial conditions, are presently 
harboring similar numbers of Class protostars. This suggests 
that after the onset of gravitational collapse, the detailed manner 
in which the collapse has been initiated does not affect much 
the ability of the clump to form stars. Moreover, both W40 
and Serpens South are characterized by large SFRs, an order 
of magnitude higher than the typical SFRs observed for embed- 
ded infrared clusters (Lada & Lada 2003). This suggests that 
(sub)millimeter protoclusters are star forming clumps caught in 
the very act of efficiently forming stars, while embedded infrared 
clusters may already be in the decay phase following the ini- 
tial burst of star formation. If they last long enough, such high 
SFRs are likely to produce a large number of stars before gas 
dispersal processes (outflows, radiation, etc.) from the growing 
cluster membership terminate the process (see the conclusions 
of Nakamura et al. 201 lb| regarding Serpens South), i.e reach 
SFEs >30% which would allow them to remain bound (Lada & 
Lada|2003|rBaumgardt & Kroupa||2007) |Goodwin||2009[ ). W40 



and Serpens South are therefore different but good candidates to 
form rich bound clusters. 
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Fig. 6. The table listing the properties of the MAMBO millimeter sources in the Aquila Rift region is available upon request to the 
first author. 
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